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A  variety  of  computational  methods,  including  the  semiempirical  techniques  AMI,  PM3,  and  MNDO,  and  the 
thermochemical  basis  sets  of  Benson  and  Stine,  was  used  to  calculate  and  compare  heats  of  formation 
(A///)  data  for  optimized  geometries  of  a  variety  of  aromatic  and  nonaromatic  hcterocyclcs.  Detailed  analyses, 
including  &-31G’  and  MP2/&-31G*  ab  initio  calculations,  were  performed  for  the  oxazolc  and  thiazolc  het- 
erocycles.  The  results  indicate  a  scatter  among  the  methods  sensitive  to  Uie  nature  of  the  heterocycle.  This 
was  in  particular  evident  in  the  oxazolc  molecule,  where  AMI  gave  a  singularly  high  value  of  A///  consistent 
with  longer  calculated  bond  lengths,  particularly  about  the  oxygen  atom.  Aromatic  stabilization  energy  appears 
to  be  addressed  differently  among  the  employed  methods.  Implications  of  this  contrast  applied  to  calculation 
of  macromolecular  systems  containing  heterocyclic  units  are  discussed. 


INTRODUCTION 

Organic  heterocycles  are  becoming  increasingly  im¬ 
portant  in  materials  chemistry.  Several  promising 
conducting  polymers,  nonlinear  optical  polymers, 
biomaterials,  and  high-strength  polymers  contain 
heterocyclic  moieties.  In  the  field  of  high-strength 
polymers,  oxazole  and  thiazole  heterocycles  are  in¬ 
corporated  into  two  advanced  high-strength  polymer 
fibers  (Fig.  1).*“'  These  materials  exhibit  extremely 
high  mechanical  strength  in  the  direction  of  the  poly¬ 
mer  backbone  due  to  the  overall  rod-like  fiber  ori¬ 
entation. 

There  is  a  great  deal  of  interest  in  determining  the 
relationship  between  computationally  derived  pa¬ 
rameters  (such  as  the  heat  of  formation)  with  ulti¬ 
mate  strength  and  elastic  modulus  properties  of 
these  polymeric  materials.  In  this  context,  we  per¬ 
formed  comparative  quantum  chemical  calculations 
on  the  cts-PBO  structure  using  the  three  most  com¬ 
monly  employed  semiempirical  methods  (MNDO,*' 
AMI,®  and  PMS).'^’®  The  RHF^  optimized  geometries 
using  MOPAC  6.0‘®  gave,  in  order  of  increasing 
A///*,  the  results  shown  in  Figure  2.  To  isolate  this 
comparison  to  the  fused  aromatic  system  (eliminat¬ 
ing  the  effect  of  interplane  angle,  «,  and  the  para- 
phenylene  unit),  these  calculations  were  repeated 
for  the  cis-benzo-bisoxazole  structure  (see  Fig.  3). 
The  magnitude  of  AA///“s  were  comparable  to  those 
seen  in  the  polymer  system,  indicating  neither  «  nor 


’Author  to  whom  all  correspondence  should  be  addressed. 


the  paraphenylene  unit  to  be  a  m:yor  contributor  to 
the  observed  differences. 

This  persistently  large  contrast  in  heat  of  forma¬ 
tion  results  led  us  to  suspect  a  significant  difference 
in  these  computational  methods  applied  to  the  ox¬ 
azole  hcterocycle.  We  thus  concentrated  our  further 
efforts  on  the  oxazole  and  the  related  thiazole  het¬ 
erocycles  to  determine  which  method  is  most  ap¬ 
propriate  for  studying  these  polymers.  (The  appli¬ 
cation  of  these  calculations  to  polymeric  systems 
will  be  addressed  in  a  subsequent  article.) 

RESULTS  AND  DISCUSSION 
Heat  of  Formation  Calculations 

Heats  of  formation  data  were  computed  for  the  ox¬ 
azole  heterocycle  using  MOPAC  6.0  and  an  eigen¬ 
vector  following  routine  for  geometry  optimization." 
The  thiazole  heterocycle,  (the  next  homologue  in  the 
chalcogen  family)  was  also  examined  due  to  the  im¬ 
portance  of  thiazole  based  high-ordered  polymers.'^ 
The  results  are  summarized  below  in  descending  or¬ 
der  of  agreement  to  experiment.  Because  experi¬ 
mental  data  for  thiazole  are  unavailable,  results  ob¬ 
tained  here  are  reported  in  the  same  order  as  that 
for  oxazole. 

A///  (kcal/mol) 


Oxazole  Thiazole 

Exp  =  -3.70  ±  0.12  Unknown 
PM3  =  -  1.48  40.28 

MNDO  =  -8.24  33.18 

AMI  =  12.46  38.63 
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tran«  poly  paraph«nylan«  beiao'bisthiazola  ftrana-PBZT^ 

Figure  1.  Structures  of  cis-poiy  paraphenylene  benzo- 
bisoxazole  (czs-PBO)  and  trans  poly  paraphenylene  benzo- 
bisthiazole  (trans  PBZT). 

AMI  predicts  a  much  higher  heat  of  formation  for 
oxazole  than  any  other  method.  This  startling  con¬ 
trast  raised  the  question  of  how  these  results  would 
compare  with  those  predicted  by  ab  initio  methods. 
This  contrast  also  supported  our  suspicion  that  the 
oxazole  heterocyclc  was  a  rntyor  source  of  discrep¬ 
ancy  between  the  scmiempirical  methods  in  the  PBO 
polymer  because  increasing  A// /  followed  the  same 
sequence  for  both  structures:  MNDO  <  PM3  <  AMI. 
For  thiazole,  the  range  of  AH /  values  is  about  half 
that  of  oxazole.  This  implies  that  the  parametric  ba¬ 
sis  sets  for  sulfur  cause  less  variation  in  the  AH/ 
results.  This  is  of  note  because  sulfur  was  only  re¬ 
cently  incorporated  into  the  AMI  model  (1990).'® 
Geometries  for  the  oxazole  and  thiazole  hetero¬ 
cycles  were  optimized  using  the  6-31G*  basis  set''' 
at  both  the  RHF  and  MP2  (Moller-Plesset  second- 
order  perturbation  theory)'®  levels  as  implemented 
in  the  Gaussian  90  program.'®*  The  resulting  molec¬ 
ular  energies  were  coupled  with  experimental  AH/s 
in  the  isodesmic  reaction  schemes  shown  in  Table  1 
to  obtain  6-31G*//6-3lG*  and  MP2/6-31G*//MP2/6- 
31G*  AH°s  for  oxazole  and  thiazole.  The  calculation 
employed  for  oxazole  and  the  ab  initio  AH/  results 
are  shown  below. 

AHj-  (oXaZOle) 

-t-  Af//(cyclopentadiene) 

-  A///(pyrrole) 

-  Ai/y’(furan) 

AHy  (oxazole)  —  AH^  j:7)(uode.sinic) 

-  A/f/(cyclopentadiene) 

+  A///(pyrrole) 

-I-  A///(furan) 

Oxazole  Thiazole 

AH/  AH/ 

Method  (kcal/mol)  Method  (kcal/mol) 

Exp  -  3.70  Exp  Unknown 

MP2/6-31G*  -4.13  MP2/6-31G*  35.54 

6-31G*  -9.25  6-31G*  32.24 


•Performed  on  CRAY  XMP  2/16. 


a  (deqr«««)  AHf*  (kca  1 /mol )  AAHf*  (kcal /mo)  ) 

KHDO  79.26  13.99  PM3-KNOO  AM1-PM3  AMl-KNDO 

PHI  3.24  37.98  23.99  39.38  63.37 

AMI  22.17  77.36 

Figure  2.  Interplane  angle  and  AH°  calculations/com¬ 
parisons  using  scmiempirical  metlrods  applied  to  cis-PBO. 

As  anticipated,  the  more  rigorous  ab  initio  treat¬ 
ment  (MP2/6-31G*)  gave  the  closer  to  AH/  value  for 
oxazole.  (For  thiazole,  such  a  comparison  is  not  pos¬ 
sible  because  an  experimental  value  for  thiazole  is 
not  known.)  For  the  oxazole  heterocycle,  the  PM3 
method  gives  significantly  better  agreement  with  ex¬ 
perimental  and  the  MP2/6-31G*  results  than  AMI. 
The  MNDO  results  were  intermediate  between  the 
two  ab  initio  values  for  both  oxazole  and  thiazole 
and  was  much  closer  to  the  experimental  value  than 
AMI  for  oxazole.  For  thiazole,  AMI  gives  a  better 
agreement  with  both  ab  initio  values  than  PM3. 
Here,  however,  the  magnitude  of  deviation  of  the 
AMI  and  PM3  methods  from  MP2/6-31G*  is  more 
comparable  than  for  oxazole. 

Molecular  Structural  Analysis 

To  understand  the  origin  of  this  range  of  results  in 
heats  of  formation,  a  detailed  structural  comparison 
was  carried  out  for  the  oxazole  and  thiazole  heter¬ 
ocycles  using  the  atom  numbering  scheme  shown  in 
Figure  4.  The  results  are  tabulated  in  Tables  II-V. 
Experimental  values  were  obtained  from  microwave 
spectral  analysis  of  the  gas-phase  molecules.  Due  to 
the  complete  plan^ity  of  these  molecules,  no  di¬ 
hedral  angle  data  is  presented.  A  compilation  of  per¬ 
cent  error  data  is  shown  in  Tables  VI  and  VII. 

Atomic  charge  densities  obtained  from  Mulliken 
population  analyses  were  also  compared  between 
the  predictive  methods  for  oxazole  and  thiazole  (Ta¬ 
bles  VIII  and  LX,  respectively).'^  For  both  com¬ 
pounds,  the  MNDO  method  gives  charges  consist¬ 
ently  much  closer  to  ab  initio  results  than  PM3  or 
AMI.  Also,  the  MNDO  method  places  greater  nega¬ 
tive  charge  density  on  the  N  heterocyclic  atom  than 
does  PM3  or  AMI.  For  oxazole,  greater  positive 
charge  density  is  observed  for  the  hydrogens  using 


6Hf*(kcal /mol  )  6AHy*(kca  1 /mol  ) 

KMDO  >11.19  FM3>MMOO  AM1>PM3  AMl-HNDO 

PH3  7.61  18.80  33.31  52.11 

ANl  40.92 

Figure  3.  Semiempirical  AH f  calculations/ comparisons 
applied  to  cts-benzo-bisoxazole. 
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Table  I.  Ab  initio  calculationii. 

Ab  initio  method  (isodesmic) 

Isodesmic  equation  1 


Optimization  energies  (au) 
(1  au  =  627.5085  kcal/mol) 


6,31G*22 

-208.80785 

-228.62521 

-244.63297 

- 192.79172 

8.37  X  lO  ’’  au 

MP2/6-31G“®-23 

-209.50418 

-229.33278 

-245.37077 

- 193.44967 

(  5.25  kcal/mol) 
16.52  X  10  ■*  au 

Exp  AH}  (kcal/mol) 

25.9 

-8.3 

-3.70 

32.1 

(10.37  kcal/mol) 

Isodesmic  equation  2 

6-31G* 

-208.80785 

Optimization  energies  (au) 

(1  au  =  627.5085  kcal/mol) 

-551.29035  -  567.28889 

o 

- 192.79172 

17.59  X  10-''  au 

MP2/6-31G’ 

-209.50418 

-551.95597 

-567.98762 

- 193.44967 

(11.04  kcal/mol) 
22.86  X  10  •’  au 

Exp  AH}  (kcal/mol) 

25.9 

27.6 

Unknown 

32.1 

(14..'J4  kcal/mol) 

Calculations  were  performed  on  a  CRAY  XMP  2/16. 


PM3  and  AMI  vs.  MNDO.  The  magnitude  of  the 
charges  cannot  be  easily  compared  between  the 
semiempirical  and  ab  initio  treatment  because  each 
uses  a  different  set  of  orbitals.  However,  the  general 
superior  agreement  of  the  MNDO  method  with  ab 
initio  methods  is  of  statistical  interest.  Differences 
in  the  signs  of  charge  density,  however,  can  reveal 
a  genuine  contrast  between  semiempirical  and  ab 
initio  treatment.  It  is  interesting  to  note  that  of  the 
three  semiempirical  methods  only  the  MNDO 
method  assigns  a  positive  charge  density  to  C2  in 
oxazole,  as  found  in  the  ab  initio  methods.  For  thia- 
zole,  the  signs  of  charge  density  are  consistent  for 
all  methods  for  each  atom  in  the  structure.  This 
could,  to  some  degree,  account  for  the  lower  devia¬ 
tion  in  AH/  results  observed  for  thiazole  vs.  oxazole. 

The  MNDO  method  was  consistently  much  closer 
to  the  ab  initio  results  for  bond  lengths  than  PM3 
or  AMI.  In  terms  of  percent  error  from  experimental 
value,  the  MNDO  method  had  a  somewhat  lower 
error  on  average.  Almost  all  bond  lengths  (with  the 
exception  of  the  C(5),  0  bond  with  MNDO,  6-3iG*, 


ii' 


\ 

i 

/ 


Cj— H 


'O, 


Figure  4.  Structures  and  atom  number  assignments  for 
oxazole  and  thiazole. 


and  MP2/6-31G*)  were  predicted  by  all  methods  to 
be  longer  to  varying  degrees  than  that  found  exper¬ 
imentally.  This  was  in  particular  so  for  the  AMI 
method,  which  predicted  much  longer  bond  lengths 
to  oxygen  than  that  observed  experimentally.  These 
longer  bond  distances  are  consistent  with  the  sig¬ 
nificantly  greater  AH /  obtained  for  the  oxazole  ring 
using  the  AMI  method.  Although  the  AMI  method 
incorporates  an  additioneil  set  of  parameters  in  the 
form  of  Gaussian  core-core  interaction  to  correct 
for  excessive  long-range  repulsion  in  the  original 
MNDO  core-core  repulsion  term,'*  the  MNDO 
method  gives  a  lower  AH°  and  shorter  bond  length 
values  for  oxazole.  Indeed,  bond  lengths  appear  to 
be  a  fairly  sensitive  function  of  relative  AH°  values. 
In  Table  VI,  signed  percent  errors  for  bond  lengths 
arc  arranged  in  descending  order,  paralleling  the  de¬ 
scending  order  for  the  AH/  values. 

All  methods,  other  than  AMI,  have  a  fairly  con¬ 
sistent  magnitude  of  signed  error,  although  the 
AH/  varies  over  a  7.78-kcal/mol  range  excluding  the 
AMI  value.  This  magnitude  of  error  (approximately 
9.5-11%  total)  appears  to  have  no  direct  influence 
on  agreement  of  AH /  with  experimental  (the  closest 
method  to  experimental,  MP2/6-31G*,  has  a  9.64% 
error).  For  AMI,  however,  a  significantly  higher 
signed  error  is  observed  with  a  much  higher  AH/ 
result.  The  longer  bond  lengths  of  the  AMI  method 
appear  to  be  a  meyor  structural  manifestation  of  the 
higher  AH /  result.  Any  possible  stabilization  due  to 
aromaticity  docs  not  appear  to  be  as  effectively  ad¬ 
dressed  in  the  AMI  method.  Bond  angle  errors  ap- 


78 

Table  II.  C 

SHAITCR  AND  WIERSCHKE 

irison  of  bond  lengths  by  method.  Angstroms  (%  deviation  from  experimental):  Oxazole. 

1‘  'c 

Cs  / 

H 

:,-H 

Method 

Bond 

MNDO 

AMI 

PM3 

6-3IG* 

MP2/6-3IG* 

Exp'^'*  Uncert. 

<0,002  A 

0.  C(2) 

1.3642 

1.4025 

1.3717 

1.3640 

1.3640 

1.3574 

(0.50) 

(3.32) 

(1.05) 

(0.49) 

(0.49) 

C(2),  N 

1J3368 

1.3276 

1.3257 

1.3369 

1.3374 

12915 

(3.51) 

(2.79) 

(2.65) 

(3.52) 

(3.56) 

N,  C(4) 

1.3996 

1.4053 

1.4160 

1.3991 

1.3998 

1.3954 

(0.30) 

(0.71) 

(1.48) 

(0.27) 

(0.31) 

C(4),  C(5) 

1.3905 

1.3873 

1.3683 

1.3906 

1.3911 

1.3525 

(2.81) 

(2.57) 

(1.17) 

(2.82) 

(2.86) 

C(5),  0 

1.3678 

1.3962 

1.3809 

1.3678 

1.3674 

1.3696 

(-0.13) 

(1.94) 

(0.83) 

(-0.13) 

(-0.16) 

C(2),  H(2) 

1.0862 

1.0902 

1.0900 

1.0858 

1.0861 

1.0750 

(1.04) 

(1.41) 

(1.40) 

(1.01) 

(1.04) 

C(4),  H(4) 

1.0811 

1.0896 

1.0897 

1.0812 

1.0814 

1.0751 

(0.56) 

(1.35) 

(1.35) 

(0.56) 

(0.59) 

C(5).  H(5) 

1.0817 

1.08.34 

1.0855 

1.0818 

1.0819 

1.0732 

(0.80) 

(0.96) 

(1.14) 

(0.80) 

(0.81) 

pear  to  be  less  systematic.  Indeed,  the  PM3  angle  : 

angles  PM3  gives  the  highest  or  lowest  predicted 

error  is  the  highest  of  all  methods  used,  although  it  ; 

angle  and  in  alternating  fashion  [for  the  C(4),  C(5), 

came  in  second  only  to  MP2/6-31G*  method  for  ac-  0  angle,  all  predictions  were  within  a  narrow  margin 

curacy  of  A///.  Consistently,  for  the  interior  ring  of  the  experimental  value].  All  methods  predicted 


Table  III.  Comparison  of  bond  angles  by  metliod,  degrees  (_%  deviation  from  experimental):  Oxazole. 
H  X, 


Method 


Angle 

MNDO 

AMI 

PM3 

621G* 

MP2/6-31G* 

Exp2< 

0,  C(2),  N 

112.40 

112.42 

110.06 

112.46 

112.43 

114.99 

(-225) 

(-223) 

(-4.29) 

(-2.20) 

(-222) 

C(2),  N,  C(4) 

10529 

105.83 

107.90 

105.27 

10520 

103.92 

(1.32) 

(1.84) 

(3.83) 

(1.30) 

(1.24) 

N,  C(4),  C(5) 

10821 

109.08 

106.48 

10822 

10829 

109.04 

(-0.76) 

(0.03) 

(-224) 

(-0.76) 

(-0.68) 

C(4),  C(5),  0 

107.73 

107.62 

108.58 

107.75 

107.65 

108.14 

(-0.38) 

(-0.48) 

(0.40) 

(-0.36) 

(-0.45) 

C(5),  0.  C(2) 

10626 

105.05 

106.97 

106.31 

106.42 

103.91 

(226) 

(1.10) 

(2.95) 

(2.31) 

(2.41) 

0,  C(2),  H(2) 

11927 

116.11 

119.01 

119.42 

119.37 

117.09 

(1.95) 

(-023) 

(1.64) 

(1.99) 

(1.95) 

N,  C(2),  H(2) 

12822 

131.46 

130.92 

128.12 

12820 

127.92 

(023) 

(2.77) 

(225) 

(0.16) 

(0.22) 

N,  C(4),  H(4) 

121.45 

120.99 

121.77 

121.49 

121.41 

121.89 

(-026) 

(-0.74) 

(-0.10) 

(-0.33) 

(-0.39) 

C(5),  C(4),  H(4) 

13024 

129.93 

131.75 

130.30 

130.29 

129.07 

(028) 

(0.67) 

(2.08) 

(0.95) 

(0.95) 

C(4),  q5),  H(6) 

134.12 

13625 

135.61 

134.12 

134.17 

135.00 

(-0.66) 

(1.45) 

(0.45) 

(-0.66) 

(-0.61) 

0,  C(5),  H(5) 

118.15 

115.42 

115.57 

118.13 

118.17 

11626 

(1.10) 

(-123) 

(-1.10) 

(1.09) 

(1.12) 

COMPARISON  OF  COMPUTATIONAL  METHODS 
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Table  IV.  Comparison  of  bond  lengths  by  method,  Angstroms  (T.  deviation  from  experimental):  Thiazolc. 


H 

c.  \ 

II  C3- 

H  S, 

Bond 

-H 

Method 

MNDO 

AMI 

PM3 

6-31G* 

MP2/6-31G* 

Exp2** 

S,  C(2) 

1.6856 

1.7109 

1.7482 

1.6860 

1.6874 

1.7239 

(-2.22) 

(-0.76) 

(1.41) 

(-2.20) 

(-2.12) 

±0.0009 

C(2),  N 

1.3240 

1.3284 

1.3236 

1.3236 

1.3232 

1.3042 

tl-52) 

(1.86) 

(1.49) 

(1.49) 

(1.45) 

±0.0011 

N,  C(4) 

1.4034 

1.3877 

1.4090 

1.4035 

1.4024 

1.3721 

(2.28) 

(1.14) 

(2.69) 

(2.29) 

(2.21) 

±0.0002 

C(4),  C(5) 

1.3788 

1,3880 

1.3679 

1.3788 

1.3795 

1.3670 

(0.87) 

(1.53) 

(0.06) 

(0.86) 

(0.91) 

±0.0004 

C(5),  S 

1.6739 

1.6761 

1.7284 

1.6746 

1.6742 

1.7130 

(-229) 

(-2.16) 

(0.90) 

(-2.24) 

(-2.27) 

±0.0003 

C(2).H(2) 

1.0861 

1.0931 

1.0909 

1.0858 

1.0860 

1.0767 

(0.87) 

(1.52) 

(1.31) 

(0.84) 

(0.86) 

±0.0018 

C(4),  H(4) 

1.0859 

1.0958 

1.0919 

1.0862 

1.0860 

1.0798 

(0.56) 

(1.48) 

(1.12) 

(0.60) 

(0.57) 

±0.0001 

C(5),  H(5) 

1.0794 

1.0869 

1.0887 

1.0795 

1.0792 

1.0765 

(027) 

(0.97) 

(1.13) 

(0.28) 

(0.25) 

±0.0002 

bond  angles  for  the  heteroatoms,  0,  and  N,  to  be 
slightly  greater  than  that  determined  experimentally. 

Critical  computational  method  analysis  for  the 
thiazole  heterocycle  is  hampered  by  the  lack  of  an 


experimentally  determined  MI/  in  the  literature. 
However,  a  close  practical  estimation  of  an  experi¬ 
mental  MI/  for  thiazole  can  be  obtained  by  com¬ 
bining  the  experimental  MI/  for  4-methyl  thiazole'** 


Table  V.  Comparison  of  bond  angles  by  method,  degrees  (%  deviation  from  experimental):  Thiazole. 


H 


\  /N3 

C.  % 


H 


■s, 


/ 


Cj-H 


Method 


Angle 

MNDO 

AMI 

PM3 

621G‘ 

MP2/6-31G* 

Exp25 

S,C(2),N 

115.07 

114.98 

114.16 

115.00 

114.98 

115.18 

(-0.10) 

(-0.18) 

(-0.77) 

(-0.15) 

(-0.17) 

±0.01 

C(2),  N,  C(4) 

109.46 

109.20 

111.55 

109.49 

109.48 

110.12 

(-0.60) 

(-0.84) 

(1.30) 

(-0.57) 

(-0.58) 

±0.02 

N,  C(4),  C(5) 

113.86 

114.92 

113.38 

113.85 

113.99 

115.81 

(-1.69) 

(-0.77) 

(-2.10) 

(-1.69) 

(-1.57) 

±0.01 

C(4),  C(5),  S 

110.39 

110.22 

111.85 

110.39 

110.26 

109.57 

(0.75) 

(0.59) 

(2.08) 

(0.75) 

(0.63) 

±0.01 

C(5),  S,  C(2) 

91.22 

90.69 

89.06 

91.27 

91.34 

89.33 

(2.12) 

(1.52) 

(-0.30) 

(2.17) 

(2.24) 

±0.03 

S,  C(2),  H(2) 

122.14 

123.23 

12426 

122.13 

122.19 

121.26 

(0.73) 

(1.62) 

(2.47) 

(0.72) 

(0.77) 

±0.5 

N,  C(2),  H(2) 

122.79 

121.79 

121.58 

122.87 

122.87 

123.56 

(-0.62) 

(-1.43) 

(-1.60) 

(-0.56) 

(-0.56) 

±0.5 

N.C(4),H(4) 

118.42 

118.66 

118.93 

118.40 

11829 

119.35 

(-0.78) 

(-0.58) 

(-025) 

(-020) 

(-0.89) 

±0.01 

C(5),  C(4),  H(4) 

127.72 

126.42 

127.69 

127.75 

128.73 

124.84 

(2.31) 

(1.26) 

(228) 

(2.33) 

(3.11) 

±0.01 

C(4),  C(5),  H(5) 

127.77 

125.58 

125.45 

127.60 

127.87 

129.03 

(-0.97) 

(-2.67) 

(-2.78) 

(-1.11) 

(-0.90) 

±0.03 

S,  C(5),  H(5) 

121.83 

124.20 

122.71 

122.01 

121.87 

121.40 

(0.36) 

(2.31) 

(1.08) 

(0.51) 

(0.39) 

±0.03 
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Table  VI.  Compilation  of  percent  error:  Oxazole. 


Method 

Bond  lengths 

Bond  angles 

AH/ 

(kcal/mol) 

Unsigned 

Average* 

Signed 

Average* 

Unsigned 

Average'’ 

Signed 

Average'’ 

AMI 

15.06 

1.88 

15.06 

1.88 

13.36 

1.21 

2.34 

0.21 

12.46 

PM3 

11.06 

1.38 

11.06 

1.38 

21.54 

1.96 

5.87 

0.53 

-1.47 

MP2/6-31G’ 

9.80 

1.23 

9.64 

1.21 

12.24 

1.11 

3.54 

0.32 

-4.13 

MNDO 

9.65 

1.21 

9.51 

1.19 

12.34 

1.12 

3.54 

0.32 

-8.24 

6-31G* 

9.60 

UO 

9.47 

1.18 

12.09 

1.10 

3.49 

0.32 

-9.25 

•Averaged  over  all  eight  bond  distances,  Table  II. 
'’Averaged  over  all  11  apgles,  Table  III. 


and  the  enthalpic  additivity  value  for  the  methyl  sub¬ 
stituent  from  Benson’s  method'®: 

Af//(thiazole)  =  AH/(4-methylthiazole) 

-  Aff/(methyl) 

=  26.7  -  (-10.08)  (Real/ mol) 

=  36.78  (kcal/mol) 

In  Table  VII,  signed  percent  errors  for  bond  lengths 
are  again  arranged  in  descending  order  and  com¬ 
pared  with  A///  values. 

The  spread  of  Afif/  values  for  thiazole  (8.04  Real/ 
mol)  is  considerably  less  than  that  of  oxazole  (21.71 
Rcal/mol).  An  examination  of  compiled  percent  er¬ 
rors  for  the  bond  length  structural  parameter  (Table 
VII)  shows  a  comparable  unsigned  percent  for  all 
methods.  PM3,  however,  has  the  greatest  signed  er¬ 
ror  of  all  methods  (identical  to  its  unsigned  error), 
indicating  that,  for  thiazole,  the  PM3  method  cal¬ 
culates  too  long  a  bond  length  of  all  bonds  in  the 
molecule.  Consistent  with  oxazole,  where  AMI  pre¬ 
dicted  the  longest  bond  lengths  and  the  highest 
A///,  so  here  does  PM3  calculate  the  highest  MI°. 
Again,  bond  angle  errors  are  less  systematic  and  do 
not  appear  to  have  any  significant  effect  on  Mf/ 
values. 

The  thiazole  vs.  oxazole  comparison  points  to 
ami’s  anomalous  behavior  for  the  oxazole  hetero¬ 
cycle  to  be  directly  related  to  the  oxygen  atom — the 
only  structural  variant  in  this  contrast.  Indeed,  if  the 
average  signed  percent  error  for  the  0,  C(2)  and  C(5), 
0  bond  lengths  is  determined  for  each  method  AMI 
gives  the  greatest  positive  error  by  a  substantial  mar¬ 


gin:  AMI  (2.63%),  PM3  (0.94%),  MNDO  (0.19%),  6- 
31G'  (0.18%),  and  MP2/6-31G‘  (0.17%).  PM3’s  treat¬ 
ment  of  the  sulfur  atom  presumably  accounts  for  its 
high  value  of  AH  °  for  thiazole,  but  here  the  contrast 
between  the  methods  is  of  much  less  magnitude  than 
for  oxazole.  The  effect  of  sulfur  parameterization 
differences  appears  to  be  much  less. 

A  structural  comparison  between  the  semiempir- 
ical  methods  focused  on  bond  length  was  also  ap¬ 
plied  to  the  larger  fused  ring  aromatic  system,  cis- 
benzo-bisoxazole.  The  data  is  presented  in  Table  X, 
Lengthening  of  the  C(l),  0(2)  and  C(1 1),  0(10)  bonds 
with  AMI  vs.  MNDO  or  PM3  was  comparable  to  that 
seen  in  the  oxazole  system.  In  the  fused  benzene 
ring  system,  there  is  a  comparable  elongation  of  the 
C — C  bonds  common  to  the  oxazole  system:  C(3), 
C(4)  and  C(7),  C(8)  for  MNDO  (0.046  A)  and  AMI 
(0.055  A).  PM3’s  bond  elongation  is  less  (0.033  A). 
However,  the  large  difference  between  the  MNDO 
and  AMI  AH  °  values  (Fig.  3)  despite  this  comparable 
phenyl  bond  elongation  indicates  this  structural  fac¬ 
tor  not  to  be  important  in  accounting  for  the  AAH  °s. 

An  identical  structural  comparison  was  performed 
for  the  cis-PBO  repeat  unit.  Completely  analogous 
bond  length  comparisons  were  seen  in  the  fused 
aromatic  portion  of  this  structure.  For  the  para- 
phenylene  system,  essentially  identical  (and  equiv¬ 
alent  around  the  ring)  bond  lengths  were  observed 
for  all  three  methods.  These  comparisons  under¬ 
score  the  prime  structural  manifestation  of  AMl’s 
high  AH/  to  be  related  to  the  oxygen  atom.  While 
the  greater  C,  0  bond  lengths  of  AMI  may  not  in 
themselves  generate  the  higher  AH/,  this  bond 


Table  VII.  Compilation  of  percent  error:  Thiazole. 


Bond  lengths 

Bond  angles 

AH/ 

(kcal/mol) 

Method 

Unsigned 

Average* 

Signed 

Average* 

Unsigned 

Average'’ 

Signed 

Average'’ 

PM3 

10.11 

1.26 

10.11 

1.26 

17.11 

1.56 

1.31 

0.12 

4028 

AMI 

11.41 

1.43 

5.59 

0.70 

13.77 

1.25 

0.84 

0.08 

38.63 

6431G‘ 

10.79 

1.35 

1.92 

0.24 

11.35 

1.03 

1.59 

0.14 

3224 

MNDO 

10.88 

1.36 

1.87 

0.23 

11.03 

1.00 

1.50 

0.14 

33.18 

MP2/6431G* 

10.64 

1.33 

1.86 

0.23 

11.82 

1.07 

2.46 

022 

35.54 

•Averaged  over  all  eight  bond  distances,  Table  IV. 
•’/.veraged  over  all  11  bond  angles.  Table  V. 
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Tlible  VIII.  Atomic  charge  densities:  Uxazole. 


Atom 

Method 

MNDO 

AMI 

PM3 

6-31G‘ 

MP-2/6-31G* 

0 

-0.1375 

-0.1369 

-0.1373 

C(2) 

0.0687 

0.0686 

N 

-0.2248 

-0.1423 

-0.1291 

-0.2246 

-0.2242 

C(4) 

-0.1745 

-0.1524 

-0.0665 

-0.0671 

C(5) 

-0.1203 

-0.0884 

-0.0484 

-0.0479 

H(2) 

0.1566 

0.2249 

0.1715 

0.1565 

0.1566 

H(4) 

0.1187 

0.1907 

0.1523 

0.1187 

0.1186 

H(5) 

0.1325  ' 

0.2031 

0.1582 

0.1325 

0.1327 

lengthening  can  be  described  as  a  discernible  struc¬ 
tural  effect  of  AMI'S  parameterization  of  oxygen 
when  applied  to  an  aromatic  system. 

From  Figures  2  and  3,  it  can  be  seen  that  AA///s 
between  semiempirical  methods  for  the  cis-PBO  re¬ 
peat  unit  and  the  cis-benzo-bisoxazole  system  are  of 
comparable  magnitude.  (The  slightly  greater  AA/// 
for  cis-PBO  can  be  related  to  minor  influences  of 
polymer  vs.  monomer  geometiy  optimizations,  in¬ 
terplane  angle  («),  and  the  paraphenylene  unit.) 
Clearly,  however,  the  fused  aromatic  structure  ac¬ 
counts  for  most  of  the  AAH  °s  observed. 

These  AA/f/s  from  Figure  3  are  considerably 
greater  than  those  observed  for  the  oxazole  structure 
(PM3-MNDO  =  7.77  kcal/mol,  AM1-PM3  =  13.94 
kcal/mol).  However,  the  PM3-MNDO  AAH/  for  the 
CIS  benzo-bisoxazole  system  can  be  largely  ac¬ 
counted  for  by  structural  constituents  of  the  fused 
aromatic  system  as  follows: 

Structural  contribution 
(PM3  vs.  MNDO) 

Two  oxazole  units  (7.77  x  2)  = 

One  benzene  ring 
(PM3-MNDO  A///)  = 

(Table  XI) 

This  value  is  close  to  the  PM3-MNDO  AAH/  ob¬ 
served  for  the  cis  benzo-bisoxazole  system  (18.80 
kcal/mol). 


The  AM1-PM3  AAH°  is  more  difficult  to  address. 
A  calculation  in  the  manner  just  shown  gives  the 
following: 

Structural  contribution  (AMI  vs.  PM3)  kcal/mol 
Two  oxazole  units  (13.94  x  2)  =  27.88 

One  benzene  ring  (AM  I-PM3  A// °)  =  -0.50 

(Table  XI)  27.33 

This  does  not  as  completely  account  for  the  AMl- 
PM3  AAH°  observed  for  cis-benzo-bisoxazole,  33.32 
kcal/mol.  The  greater  A  AH/  discrepancy  could  be 
due  to  the  compounding  effect  of  AMl’s  relative  min¬ 
imization  of  aromatic  stabilization  energy  applied  to 
a  larger  aromatic  structure. 

Comparison  with  Thermochemical  Data  Bases 

The  oxazole  heterocycle  appears  to  be  a  .structure 
maximizing  the  predictive  differences  of  semiempir¬ 
ical  and  ab  initio  methods.  This  raises  a  question  of 
how  these  semiempirical  methods  contrast  with 
each  other  and  experimental  heats  of  formation  for 
a  range  of  heterocyclic  compounds.  The  question  of 
aromatic  stabilization  of  oxazole  obviously  being 
addressed  very  differently  by  the  AMI  method 
prompted  a  compilation  of  results  for  two  classes 
of  heterocyclic  compounds:  aromatic  and  nonaro¬ 
matic. 


kcal/mol 

15.54 

1.20 

16.74 


Table  IX.  Atomic  charge  densities;  Thiazole. 

Method 

Atom  MNDO  AMI  PM3  6-31G*  MP-2/6-31G‘ 


Table  X.  Comparison  of  bond  lengths  for  semiempirical  methods;  cis-Benzo-Bisoxazole. 
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'Obtained  from  ref.  26. 

‘'Obtained  from  ref.  21. 

'Derived  from  Aff/  for  4-methylthiazole  (ref.  18).  This  compound  is  therefore  not  included  in  comparative  study. 
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Atoa  1 

Notat ion 

Coaputat ion 

Total 

(kcal/mol ) 

(kcal /mo 

1.11 

C(M2,C1.01) 

2 

X  -2.96 

-5.92 

2,10 

0(C1,C1) 

2 

X  12.42 

24.85 

5,12 

N(C2,C1) 

2 

X-16.94 

-33.88 

4.8 

C(H1,C5,C5) 

2 

X  20.74 

41.48 

3,7 

C(Ol,C5.C5) 

2 

X-16.28 

-32.56 

6,9.14,15,17,18 

C(C5,C5,H1) 

6 

X  3.27 

19.62 

13,16 

C(C1,C5,C5) 

2 

X  5.41 

10.82 

Strain  corractions 

1,3.4,11 

C<1S(2.1) 

4 

X  0.62 

2.46 

2,10 

05(1.1) 

2 

X  0.95 

1.90 

5,12 

N5(2. 1) 

2 

X  -3.85 

-7.70 

7,8 

C5{l.l) 

2 

X  1.19 

23.45 

*• — 

'53  (  1  14 

\W/ 

14  — «S 


Atoa  I  NotAtion 

1. 11 
2,10* 

5,12 
4,B 
3,2 
13,14 

6,9, 14.15, 17, li 
Ring  Strain  Corr«ction 


Cd-(O) (C; 

O-(Cd) (Cb) 

»I-(Cb) 

Cb-(H) 

Cb-<0) 

CB-(Cd) 

Cb-(M) 


Co«put  at  1  on 
(kcal /mo\ ) 
10.3  M  2 
-27.05  X  2 
16.7  V  2 
-0.5  K  2 
-0.9  X  2 
5.68  X  2 
3.30  X  6 


Tot  a  1 

(kcal /aol ) 
20.6 
-54 . 1 
33.4 
-1.0 
-1.8 
11.36 
19.8 


-5.8  X  2 


-11.6 

16.66 


*K«kul«  stnictur*  usad  to  assign  strain  corractions  for  fused 
aromatic  structura. 


*  0-(Cd){CB) 

-21.1 


is  an  average  of  0-(C4)2 


-33.0  and  0-<Cb)2  • 


Figure  5.  Stine/ Kramer  MJ/  calculation  for  cis-PBO. 

Two  thermochemical  data  bases  developed  by 
Benson  and  Stine/Kramer  were  also  employed  to 
determine  heats  of  formation  in  this  comparative 
study.^  The  recently  developed  Stine/ Kramer 
method  (1989)  employs  constituent  A///  values  for 
all  possible  atomic  bonding  schemes  of  compounds 
containing  any  combination  of  C,  H,  N,  and/or  O. 
These  values  were  obtained  empirically  by  least- 
squares  analysis  of  gas-phase  data  of  1 129  molecules. 
This  approach  differs  from  earlier  thermochemical 
data  bases  (i.e.,  Benson’s)*'*  in  that  it  ( 1)  is  structured 
toward  energetic  molecules  and  (2)  uses  fewer  pa¬ 
rameters  that  span  a  larger  range  of  organic  mole¬ 
cules. 

Calculations  using  the  Stine/ Kramer  and  Benson 
methods  to  determine  AH/  values  involve  summa¬ 
tion  of  constituent  values  for  specific  bonding  en¬ 
vironments  of  each  atom  in  the  molecule  of  interest, 
plus  strain/ interaction  energy  correction  terms.  A 
different  notation  scheme  is  used  in  the  Stine /Kra¬ 
mer  method  to  avoid  confusion. 

A  detailed  accounting  for  the  Stine/Kramer  AH/ 
calculation  for  the  cis-PBO  and  oxazole  structures 
using  the  notation  of  ref.  20  is  presented  in  Figures 
5  and  6,  respectively.  The  uncertainty  associated 
with  the  Stine/Kramer  method  is  ±5-6  kcal/mol.*“ 

For  tliC  Benson  calculations,  several  best  matches 
had  to  be  derived  from  similar  bonding  environments 


**  Closaat  valua  froa  table  to  parallel  oxazole  ring  systea. 

Beth  ot  thaae  are  saurcea  of  Ttor  oi  approxiaateiy  •  5 
kcal/sol . 

Figure  7.  Benson  AH/  calculation  for  cis-PBO. 

listed  in  the  tables.  This  procedure  introduces  more 
error  than  that  inherent  in  the  tabular  data  itself, 
which  overall  is  approximately  ±1  kcal/mol.'**  Ben¬ 
son  calculations  for  the  cis-PBO  and  oxazole  struc¬ 
tures  are  presented  in  Figures  7  and  8,  respectively. 

In  Table  XI,  the  proposed  experimenlcU  AH/  for 
thiazole  is  based  on  that  for  4-mcthylthiazoIe  with  a 
methyl  group  correction  from  the  Benson  method. 
However,  the  full  Benson  calculation  for  thiazole  is 
shown  in  Figure  9  and  is  used  for  the  Benson  entry 
in  Table  XI. 

The  Stine /Kramer  method  gives  a  AH°  of  23.45 
kcal/ mol  for  the  cis-PBO  system,  roughly  interme¬ 
diate  between  that  of  MNDO  and  PM3  from  Figure 
2,  while  Benson,  AH/  =  16.66  kcal/mol,  agrees 
closely  with  MNDO.  The  AM  1  result  is  still  singularly 
high. 

Further  calculations  on  the  structurally  related 
monomers  to  the  cis-PBO  system  gave  the  results 
shown  in  Figure?  10  (in  order  of  increasing  AHf°). 
Again,  the  Benson  method  agrees  with  MNDO.  The 
Stine/Kramer  method  is  intermediate  between  Ben¬ 
son  and  PM3,  and  AMI  is  considerably  higher.  For 
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Atom  f 

Notation 

Compulac ion 

Total 

kcal/mol 

(kcal /mol ) 

(kcal /mol ) 

Atom  1 

Notation 

1,3 

C<i-(0)  (H) 

8.6X2 

17.2 

1 

C(N2,01,H1) 

0.34 

2 

0-(Cd)2 

-33.0 

-3J.0 

2 

0(C1,CI) 

12.42 

4 

C<i-(Ni)  (H) 

6.78 

6.76 

3 

C(C2,01,H1} 

-17.27 

5 

Ni-(C) 

21.3 

21.3 

4 

C(M1,C2,H1) 

18.92 

S 

N(C2,C1) 

-16.94 

Ring  Strain 

Correction 

Strain  Corractions 

(!) 

1,3,4 

2 

3  X  Cd5<2,l) 
05(1,1) 

1.85 

0.95 

5 

N5(2,l) 

-3.85 

-5.8 

-3.58 

6.48 

Figure  6.  Stine/Kramer  AH/  calculation  for  oxazole. 


Figure  8.  Benson  AH/  calculation  for  oxazole. 


CUMPAl{lSt)N  OK  Cl)Ml*UTA'noNAL  MK  I  IK -OS 


Atoa  ff  Notation  Coapoi ai i on  Total 

(kcal/mol )  (kcal /mol ) 

1.3  Cd-(S)(H)  8.S6  X  2  17-12 

2  S-(Cjj)2  -4.54  -4.54 

4  C  (Hi) (H)  6.78  6.78 

4  .4  -tC)  21.3  21.3 

Ring  Strain  .  .action 


42 . 36 

Figure  9.  Benson  a///  calculation  for  thiazole. 

the  simple  oxazole  heterocyclc,  a  large  contrast  is 
again  noted  using  these  methods.  The  Stine/ Kramer 
value  of  A//^°  =  -3.58  kcal/mol  provides  the  best 
agreement  with  the  experimental  value — actually 
witliin  the  experimental  error.  This  allows  us  to  give 
the  final  ordering  of  the  oxivzole  A//,°  accuracy  for 
the  seven  methods  shown  below: 

Experimental  =  -3.70  ±  0.12  kcal/mol 
Stine/ Kramer  =  -3.58  kcal/mol 
MP2/6-31G*  =  -4.13  kcal/mol 
PM3  =  -1.48  kcal/mol 
MNDO  =  -8.24  kcal/mol 
6-31G*  =  -9.25  kcal/mol 
Benson  =  6.48  kcal/mol 
AMI  =  12.46  kcal/mol 

Here,  Benson  contrasts  sharply  with  MNDO.  For  this 
heterocyclic  compound,  PM3  is  the  best  of  the  semi- 
empirical  methods  based  upon  its  agreement  with 
the  experimental  value. 

To  assess  the  general  utility  of  these  techniques 
in  treating  heterocyclic  compounds,  the  MNDO, 
PM3,  AMI,  Stine/Kramer,  and  Benson  methods  were 
used  to  calculate  AH/  for  optimized  geometries  of 
a  series  of  aromatic  and  nonaromatic  heterocycies. 
(The  Stine/Kramer  and  Benson  geometries  arc  in¬ 
herent  in  the  additivity  tables  employed.)  For  the 
semiempirical  methods,  compiled  data  were  drawn 
fi’om  the  article  by  Stewart.**  Compounds  not  in¬ 
cluded  in  that  compilation  were  calculated  using 


MHOO  -  11.45  kcal/BOl  KMDO  *  3.64  kcal/nol 
B«nson  >  12.58  kc«l/«ol  Banson  •  4.58  kcal/»ol 
Stina  •  21.00  kcal/aol  Stina  -  19.28  kcal/mol 
PH3  •  33.47  kcal/mol  PH3  •  23.90  kcal/mol 
ANl  -  68.00  kcal/mol  AMI  •  60.27  kcal/mol 

Figure  10.  AH/  calculations  for  monomer  systems. 


8.'j 

MOPAC  6.0.*"  (The  Stine/ Kramer  method  is  not  ap¬ 
plicable  to  compounds  containing  sulfur.)  These  re¬ 
sults  were  compared  with  experimental  Allf"  values 
in  the  gaseous  phase-'  and  are  shown  on  Tables  .\1 
and  Xll  for  aromatic  and  nonaromatic  compounds, 
respectively.  Because  a  true  experimental  A//,°  for 
thiazole  is  not  known,  this  compound  was  excluded 
from  this  final  comparative  .study.  (In  Table  XI,  ben¬ 
zene  is  included  as  a  contrast  with  homocyclic  .sys¬ 
tems.) 

To  maintain  as  much  simplicity  as  im.ssible  in 
these  comparisons,  structural  features  beyond  tliat 
of  the  heterocycle  unit  (i.e.,  functional  groups  or 
alkyl  appendages)  have  in  general  not  been  included. 
Exceptions  to  this  are  the  dibenzo  derivatives  of 
pyrrole,  furan,  and  thiophene,  the  benzo  derivative 
of  oxadiazolc,  and  methyl  substitution  of  the  thiazole 
and  furazan  rings.  The  dibenzo  derivatives  were  used 
because  of  their  availability  for  all  three  N,  O,  and 
S  heterocycies  and  so  provided  a  convenient  com¬ 
parison  to  a  larger  structural  system.  The  other  mod¬ 
ified  structures  were  used  because  these  were  the 
simplest  compounds  of  their  class  with  reported 
AH°  data  in  the  gaseous  state. 

The  degree  of  agreement  between  the  five  meth¬ 
ods  and  experimental  values  shows  considerable  di¬ 
versity.  However,  several  general  trends  or  patterns 
are  evident: 

1.  For  most  aromatic  compounds  (with  the  excep¬ 
tion  of  thiophene,  dibenzothiophene,  furan,  and 
pyridine),  the  Stine/ Kramer  or  P.V13  methods 
place  first  or  second  in  agreement  with  experi¬ 
mental  values.  The  Bemson  and  Stine/ Kramer 
methods  collectively  work  best  for  the  simple 
monoheteroatomic  saturated  heterocycies.  P.M3, 
overall,  is  the  best  semiempirical  method  for  this 
class  of  compounds. 

2.  For  most  of  the  N-containing  six-membered  aro¬ 
matic  ring  compounds,  AM  1  gives  the  best  semi¬ 
empirical  result.  However,  for  pyridazine  (adja¬ 
cent  N  atoms)  and  for  most  five-membered 
aromatic  heterocycies  of  one  nitrogen  and/or  one 
oxygen  AM  1  gives  the  poorest  agreement,  usually 
by  a  substantial  positive  error. 

3.  While  MNDO  docs  most  poorly  for  the  N-contain¬ 
ing  six-membered  aromatic  compounds,  it  is  the 
best  semiempirical  method  for  several  key  five- 
membered  aromatic  hcterocyclics;  furan,  thio¬ 
phene,  isoxazole,  and  their  dibenzo  derivatives. 

4.  PM3  gives  the  best  semiempirical  results  for  ox¬ 
azole,  4-methylthiazole,  and  pyrrole. 

5.  For  compounds  having  three  or  more  hetero¬ 
atoms  in  connectivity,  there  is  a  wide  scatter 
among  the  methods. 

6.  The  Stine/ Kramer  metliod  compares  quite  favor¬ 
ably  with  the  three  semiempirical  methods  in  its 
agreement  with  experimental  AH°  values.  For  ar¬ 
omatic  compounds,  the  Stine/ Kramer  method 
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Table  XIII.  Frequency  of  closesl  agreement  with  experimental  A///’. 


Method 

Compound  class 

Stine/Kramer 

Benson 

MNDO 

AMI 

PM3 

Aromatics 

e 

5 

2(8) 

0(3) 

4(6) 

Overall  average  error 
(unsigned;  from  Table  IX) 

3.34 

7.50 

8.68 

13.21 

7.34 

Nonaromatics 

2 

8 

0(2) 

0(3) 

2(7) 

Overall  average  error 
(unsigned;  from  Table  X) 

3.28 

1.G4 

10.85 

8.55 

5.08 

gives  the  greatest  agreement.  For  nonaromatics, 

Benson  gives  the  best  results. 

Table  Xlli  summarizes  these  results  giving  the 
number  of  times  each  method’s  computed  MI/  was 
closest  to  the  experimental  value  and  its  overall  un¬ 
signed  error.  The  same  comparison  for  only  the  semi- 
empirical  techniques  (MNDO,  AMI,  PM3)  appears  in 
parentheses. 

In  comparing  average  errors  for  all  methods, 
clearly  Stine/Kramer  is  superior  for  aromatic  com¬ 
pounds.  It  shows  the  lowest  overall  error  in  terms 
of  absolute  magnitude.  For  the  aromatic  compounds, 
Stine/Kramer  is  followed  by  PM3,  Benson,  MNDO, 
and  AMI  in  terms  of  average  unsigned  error.  For  the 
nonaromatics,  Benson  is  followed  by  Stine/ Kramer, 
PM3,  AMI,  and  MNDO. 

Among  the  semiempirical  techniques,  it  appears 
the  PM3  and  MNDO  methods  best  address  aromatic 
systems.  AMI  is  the  best  method  only  for  specific 
N-containing  six-membered  ring  compounds.  This  is 
a  surprising  result  because  the  MNDO  method  is  the 
oldest  of  these  semiempirical  methods  and  has  the 
most  limited  basis  set  of  empirical  parameters.  How¬ 
ever,  for  nonaromatic  compounds  a  dramatic  rever¬ 
sal  of  MNDO’s  performance  is  noted.  Here,  it  gives 
the  worst  agreement  with  experiment,  consistently 
with  negative  error.  For  this  compound  class,  PM3 
gives  best  agreement  of  the  three  semiempirical 
methods.  For  both  classes  of  compounds,  AMI  gives 
the  poorest  overall  agreement  with  experimental 
data. 

From  this  comparison  of  aromatic  vs.  nonaromatic 
heterocyclics,  it  appears  that  the  MNDO  parameter¬ 
ization  incorporates  as  effectively  as  PM3  the  inher¬ 
ent  aromatic  stabilization  in  the  predicted  MI /  val¬ 
ues.  With  the  exception  of  the  six-membered 
N-containing  aromatic  heterocyclics,  where  AMI 
gave  an  average  error  of  -4.25  kcal  (vs.  -7.83  for 
PM3,  - 12.45  for  MNDO,  and  - 1.70  for  Stine/Kra¬ 
mer),  the  AMI  method  consistently  gave  a  large  pos¬ 
itive  error  for  the  MI/  value  for  aromatic  com¬ 
pounds.  Aromatic  stabilization  energy  appears  to  be 
less  adequately  addressed  by  the  AMI  method  for 
five-membered  heterocyclics  in  general.  The  con¬ 
sistently  longer  bond  lengths  in  the  AMl-optimized 
geometry  appears  to  underscore  the  excessive  im¬ 
portance  AMI  places  on  repulsive  atomic  interac¬ 


tions.  This  factor  would  tend  to  negate  any  stabili¬ 
zation  (lowering  of  bond  lengths  and  MI /)  due  to 
aromaticity. 


CONCLUSION 

We  conclude  that  the  Stine/Kramcr  and  Benson 
methods  give  the  best  overall  agreement  with  ex¬ 
perimental  MI/s  where  applied  to  aromatic  and 
nonaromatic  heterocyclic  compounds,  respectively. 
The  empirically  derived  constituent  basis  set  of  the 
Stine/Kramer  method  focuses  on  applicability  to  en¬ 
ergetic  compounds.  Approximately  30%  of  energetic 
compounds  could  be  classified  as  aromatic.  (Unfor¬ 
tunately,  the  Stine/ Kramer  model  is  limited  by  its 
not  considering  hypervalent  compounds  and  only 
those  containing  C,  H,  N,  and  0.)  Ring  strain  energy, 
represented  in  the  study  here  by  the  three-  and  four- 
membered  nonaromatic  heterocyclics,  is  best  ad¬ 
dressed  by  thf  Benson  method. 

The  PM3  method,  parameterized  from  763  com¬ 
pounds  (106  of  which  are  hypervalent),  comes  clos¬ 
est  to  the  Stine/Krauner  and  Benson  methods  of  the 
three  semiempirical  models.  Its  scope  is  broader 
than  Stine/Kramer  and  for  Utis  reason  it  provides  an 
effective  complement  to  the  Stine/Kramer  method 
for  heterocyclic  compounds.  It  is  remarkable  how 
well  MNDO  does  for  aromatic  heterocyclics,  consid¬ 
ering  only  34  molecules  were  used  to  develop  the 
C — H — N — 0  set  and  up  to  a  few  tens  of  molecules 
for  the  other  elements. 

The  Stine/ Kramer  and  Benson  methods  are  not 
applicable  to  the  nonequilibrium  energy  calcu¬ 
lations  used  for  derivation  of  polymer  elastic  force 
constants  or  moduli  values.  This  study  implies  that 
each  semiempirical  method  is  best  suited  to  specific 
polymer  structures  assuming  two  points  are  valid: 
(1)  The  heterocyclic  structural  component  mostly 
accounts  for  variation  in  calculated  polymer  prop¬ 
erties  and  (2)  the  relative  accuracies  of  semiempir¬ 
ical  methods  at  highly  strained,  nonequilibrium  ge¬ 
ometries  parallels  that  seen  for  the  ground-state 
equilibrium  geometry.  From  this  perspective,  the  cis- 
PBO  structure  and  others  like  it  having  fused-ring 
five-membered  heterocyclics  are  best  treated  overall 
by  the  PM3  method.  On  the  other  hand,  N-containing 
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six-membered  heterocyclics  (for  example,  the  qui- 
noxaline  polymers)  are  best  treated  by  AMI.  Even 
MNPO  may  provide  a  comparably  accurate  treat¬ 
ment  of  furan,  thiophene,  and  isoxozole-based  poly-- 
mers.  Overall,  PM3  performs  best  of  the  semiempir- 
ical  methods,  but  there  is  a  'vide  scatter  sensitive  to 
the  nature  of  the  heterocycle.  Such  a  diversity  of 
results  obtained  by  the  methods  employed  here  for 
a  fundamental  property  such  as  heat  of  formation 
implies  the  danger  of  an  even  more  exaggerated  scat¬ 
ter  for  predicted  properties  of  a  macromolecular  ar¬ 
ray  (i.e.,  polymers).  Of  particular  concern  in  the  poly¬ 
mer  area  is  the  effect  of  this  method  diversity  on 
prediction  of  ultimate  physical  property  terms  such 
as  elastic  moduli.  An  extremely  useful  future  pre¬ 
dictive  tool  for  a  polymer/ material  chemist  would 
involve  a  table  of  force  constants  for  specific  poly¬ 
mer  structural  units.  This  table  could  then  be  applied 
to  any  conceivable  polymeric  structure  to  predict  its 
ultimate  strength  properties.  However,  for  this  ca¬ 
pability  to  be  realized  inconsistencies  in  the  semi- 
empirical  methods  applied  to  heterocyclic  com¬ 
pounds  need  to  be  resolved.  Considering  the  central 
importance  of  aromatic  heterocyclic  compounds  in 
burgeoning  new  areas  of  materials  chemistry  (for 
example,  conducting  polymers,  NLO  materials,  and 
high-strength  fibers),  a  more  rigorous  inclusion  of 
this  class  of  compounds  into  parameterization  of 
semiempirical  methods  is  obviously  needed.  From 
this,  a  single  semiempirical  method  best  suited  to  all 
aromatic  heterocyclic  systems  may  become  avail¬ 
able.  This  would  provide  a  unified  basis  for  predic¬ 
tion  of  the  physical  and  chemical  properties  of  this 
important  class  of  compounds.  It  is  important  that 
the  implications  of  this  semiempirical  method  con¬ 
trast,  focused  here  in  terms  of  the  heterocyclic  unit, 
be  applied  to  ultimate  polymer  properties  such  as 
elastic  moduli.  These  studies  are  currently  being 
done  and  will  be  reported  in  a  subsequent  article. 
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